Exact asymptotics of the optimal L^^^^-error of linear spline 

interpolation 

Vladislav Babenko, Yuliya Babenko* Dmytro Skorokhodov 

^ Abstract 

^ In this paper we provide the exact asymptotics of the optimal weighted Lp-error, < p < oo, of 

hnear sphne interpolation of functions with positive Hessian. The full description of the behavior 
of the optimal error leads to the algorithm for construction of an asymptotically optimal sequence 

^ of triangulations. In addition, we compute the minimum of the Lp-error of linear interpolation of 
the function + over all triangles of unit area for all < p < oo. This provides the exact 

^ constant in the asymptotics of the optimal error. 

^ Keywords: spline, interpolation, adaptive, exact asymptotics, optimal error. 

g 1 Definitions, history, and main results. 
1.1 Definitions. 

^ Let be the space of points in the plane endowed with the usual Euclidian distance. The distance 
Q between points A and B in we shall denote by \AB\. 

Let D = [0, 1]^ C M?. We use this region for simplicity; the approach presented in this paper 
can be applied to any bounded connected region which is a finite union of triangles. By C{D) we 
^ shall denote the space of functions continuous on D. 

^ Let Lp[D), < p < oo, be the space of measurable functions /:£)—>■ R for which the value 

'% 11/11. =ii/iim.)-{ (/i^mi^''^) • « o<f<-. 

I, esssup{|/(T)| : r e D}, if p = oo. 

is finite. 

Remark: Note that in the case 1 < p < oo the functional || • ||p is the usual norm in the space 
Lp{D). For p e (0, 1) the above relation defines a seminorm satisfying 

(1) ll/ + ^ll^<ll/ll^ + W 

for arbitrary functions f,g E Lp{D). 
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Given a positive continuous weight function Q e C{D), we define the weighted Lp-norm || • n 

as 

ii/iu = \\m\p- 

In the case Q = 1, instead of L^ q we shall simply write Lp. 

A collection A^v = A]sf{D) = {Ti}^^ of N triangles in the plane is called a triangulation of a 
set D provided that 

1. any pair of triangles from Ajv intersect at most at a common vertex or along a common edge, 

N 

2. D = [jTi. 

1=1 

Let Pi be the set of bivariate linear polynomials p{x, y) = ax + by + c. Given a triangulation 
Ajv, define the class of linear sphnes on A to be 

S,{Aj,) := {/ e C{D) : = 1, ...,N 3p e P, such that /|t, = p|tJ. 

Let Atv) denote the spline from 5*1 (Ajv) which interpolates the function / G C{D) at the 
vertices of the triangulation A^v- Note that the linear spline Ajv) is uniquely determined by 
its values at the vertices of the triangulation A^. 

Now let the function / e C'^{D) and the number of triangles e N be fixed. Set 

-Rjv(/, Lp^^) := inf 11/ - Ajv)||p,t^, 

where inf is taken over all triangulations of D with triangles. This quantity will be called optimal 
Lp,n-error of the interpolation of the function / by splines s{f, Ajv). 
A triangulation A^^ is called optimal for the given function / if 

\\f-s{f,A%)\\p,n = RNif,Lp^a). 

Exact values of Ri\i{f, Lp^o,) as well as optimal triangulations A^ for every particular function 
/ can be found only in few exceptional situations. That is why the following two problems are 
interesting and important: 

1) Find exact asymptotics of the optimal error i?Ar(/, Lp^n) as A^" — >■ oo for any function / e 
C\D). 

2) Find an asymptotically optimal sequence of triangulations, i.e. a sequence of triangulations 
{A^}?^=i of D such that 

\\f-s{f,A*^)\\p,a ^ ^ 
N-^oo RN{f,Lp^a) 

For / e C'^{D) set /^^ := 0, f^y := fyy := 0, and denote the Hessian of / by 

H{f;x,y) := {fccxfyy - fly){x,y). 

The main purpose of this paper is to solve the two formulated above problems for functions 
/ G C^(-D) such that H{f\ x,y) > on D. We shall need to resolve the following auxiliary extremal 
problem, which is of independent interest. 
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For a given triangle T <Z D and a function g € C{D) we shall denote by l{g,T) the linear 
polynomial which interpolates g at the vertices of the triangle T. In addition, for < p < oo and 
continuous weight function fl set 



dig,T,L,,n) h - l{g,T)\\L,^,^T). 



For Q{x, y) := x"^ + y"^ define 
(2) 



C+ := inf 



d{Q, T, Lp 



where |T| denotes the area of the triangle T. 

Main results of this paper are in the following two theorems. 

Theorem 1 Let f G C'^{D); H{f;x,y) > > for all {x,y) G D. Let a positive continuous 
weight function Q{x,y) also be given. Then for all < p < oo 



lim N ■ R^if, Lp,n) 



Ci 



L p (D) 

p+i 



Theorem 2 For any p > infimum in the definition ofCp is achieved only on equilateral triangles. 
Consequently, 



\ Jil-x^-y^dxdy, 

To 



where Tq is the equilateral triangle with the center of circumscribed circle at the origin. 
Remcirk. It is easy to see that 



3^3 



1+^ 



71 f 1 

6 / x(l — x'^Y arccos — dx 

p+1 J ^ ^ 2x 



1/2 



Moreover, if B{a, b) is the Euler Beta function, and 

X 



is the incomplete Beta function, then 



(— 
V3V3 



1+^ 



TT 



P+1 



-3B B 



3 1 



2' 2 



CoroIIeiry. Let Q{x, y) — Ax'^ + By"^ + 2Cxy. Then for every p G (0, oo) 
(3) inf "^^^'^'f^^ = C;VAB - C2. 



IT 
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The idea of constructing an asymptotically optimal sequence of triangulations is to substitute 
the function / by piecewise quadratic function for every N, for which the good triangulation is 
constructed with the help of triangles solving problem ([3]). 

In the proof of the lower estimate in Theorem 1 the fact that for an arbitrary / G C'^{D) with 
positive Hessian there exists a constant k, depending on the function / only, such that 



d{f,T,L,)> .\Tr-., 



K ■ diamT 1^7^11+1 

25p+l 

played an important role (see Lemmas 7, 8, and 10). 
1.2 History. 

The first result related to piecewise linear interpolation in two dimensional case was obtained by L. 
Fejes Toth. He indicated ([7], Ch. 5, §12) that for a body C C M'^ with boundary of differentiability 
class and positive Gaussian curvature K{x, y) the Hausdorff distance of C to its best inscribed 
polytope with at most n vertices is 

1 + 0(1) f j^(^^^yY/-i(ia{x,y) 



3-\/3 \Jdc 7 

as n — )■ 00, where a is the surface area measure on dC. He also indicated that the distance of C to 
its best inscribed polytope with at most n vertices (measured as a volume of the difference between 
C and the polytope) is 



4\/3 \Jdc ' ' y n 

as n — 7- 00. Even though all the ideas were mostly contained in [7j, formally the complete proof 
was given by Gruber (see [8]). In addition, Gruber generalized these results to higher dimensions, 
however the constants were implicit. He also proved similar estimates for the error measured in 
symmetric difference metric, Banach-Mazur metric, as well as using the Schneider distance (see [^). 
Among other interesting results on these and closely related questions are results by Boroczky and 
Ludwig [31 Hj. Survey of further results on approximation of convex bodies by various polytopes 
in different metrics (inscribed, circumscribed, of the best approximation, with restrictions on the 
number of faces, etc.) can be found, for example, in Pi [S]. 

With regard to the asymptotically optimal approximation of functions by linear splines in dif- 
ferent metrics the following results are known. 

Nadler in [9] studied the sequence of asymptotically optimal triangulations for the (in general 
discontinuous) piecewise linear approximation for an arbitrary function in the sense of minimizing 
error in the L2-norm. 

In the paper [Tj this problem was solved for p = 00 and = 1, and in [2J it was extended for 
p = 00 and positive weight functions Q G C{D). 

Note that the case p = 00, i7 = 1 is close to (but not indentical with) the result of L. Fejes Toth 
on approximation of convex bodies by polygons in Hausdorff metric. The case p = 1 follows from 
the results of Boroczky and Ludwig (see |4j). 

As for the question of computing the constant and the optimality of the equilateral triangle, 
the following cases have been investigated: 
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• p = oo (D'Azevedo and Simpson, in [U]), 

• p = 1 (Boroczky and Ludwig, in [1]), 

• p = 2 (Pottmann, Hamann et al, in [lO]), 

• p = 3,4,5 (Yu. Babenko, in 0), 

• p G N (Chen, in [5]). 

The rest of the paper is organized as follows. Section 2 provides certain preliminary results, in 
particular, on how affine transformations affect the error of interpolation of a quadratic function by 
linear splines. Section 3 contains the computation of the constant in the case 1 < p < oo. The 
case p G (0, 1) is rather technical, and thus is presented in the Appendix. In Section 4 we provide 
the proof of the upper estimate in Theorem 1. This proof leads to the algorithm for construction 
of an asymptotically optimal sequence of triangulations. Section 5 contains the proof of the lower 
estimate in Theorem 1. 

2 Preliminaries. 

In order to investigate the asymptotic behavior of the optimal error of piecewise linear interpolation 
of an arbitrary function from the class C^(-D) we shall use linear interpolation of the piecewise 
quadratic functions which appear as an intermediate approximations of /. 

Some of the facts we shall present in this section are quite easy to see. However, we shall 
prove some of them, first of all for completeness, and secondly, because we shall use them in the 
construction of an asymptotically optimal sequence of triangulations. 

Let us define the modulus of continuity of g E C{D) by 

uj{g,5) ■.= ?,nY>{\g{x,y) - g{x\y')\: \x-x'\<5] \y-y'\<5] {x,y),{x',y') e D}. 
For the function / G C^{D) set 

(4) u{6) := max{uj{U^, 6), co{fyy, 6), uj{f^y, 5)}. 

Lemma 1 Let f G C'^{D). If P2 = P2{f] x,y; XQ,yQ) denotes the second degree Taylor polynomial 
of f at the point {xo,yo) inside the square C D with side length equal to h, then we have the 
following estimate: 

\\f-P2\\L^iD,)<2euj{h), 

where u}{t) is defined in Q). 

This simple lemma can be proved similarly to Lemma 1 from [1]. The following statement is 
almost obvious. 

Lemma 2 For the given quadratic function 

(5) Q{x, y) = Ax^ + By^ + 2Cxy, 
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an arbitrary triangle T, and any c G M? , the Lp-errors (0 < p < oo) of linear interpolation of 
Q{x, y) on T, c + T, and a triangle T which is symmetric to T with respect to the midpoint of any 
side ofT, are equal, i.e. 

d{Q, T, Lp) = d{Q, c + T, Lp) = d{Q, f , L^). 

For an arbitrary linear transformation S* : — ?► denote by det S the determinant of the 
matrix of this transformation. 

Lemma 3 Consider a non-singular affine mapping F = c + F , where c G and F is a linear 

rr _ |T| 

transformation. Then for any quadratic function (5) and any triangle T we have \F^'^{T) \ = ^ 

U det F 

and 



d{Q o F, F~\T), Lp) = d{Q, T, L^) ■ ( ^) ' • 

Vdet F / 



This lemma can be easily verified by a routine change of variables. 

Assume now that the function (|5| has a positive Hessian, i.e. AB — > 0. Let us find the 
eigenvalues and unit eigenvectors of the matrix of quadratic form (|5]). 
For eigenvalues we have 



A + B (A + bV .^^ ^ A + B (A + bV .^^ 

Observe that < Amin < Amax- In addition, note that 

\ \ A ID r^2 

'^min^max .rllJ O . 

Lst ^2) € be an eigenvector of Q{x, y) corresponding to the eigenvalue Amax- Then (—^2, ^1) ^ 
§^ is an eigenvector corresponding to the eigenvalue Amin- 



Lemma 4 For the quadratic form such that AB — > it follows that 

(6) d{Q,T,Lp) > C+|T|^+VAmi„Amax. 

Proof. Recall that Q{x,y) = x"^ + y"^. Obviously, Q is the canonical form of Q. We obtain this 
canonical form by the following two linear transformations Fi and F2 

(7) Fi: x = ^,x'-^2y', y = ^2x' + ^iy', 

and 

(8) F,: x' = -^, y 



Set F := Fi o F2. Note that 



detF 



V Amax mm 
1 



V^Amin An 
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and 

QoF{u,v) = u^ + w^ 
Thus, in view of Lemma [3] it follows that 

d{Q,T, Lp) =diQoF,F'' (T) , Lp) ■ (det F) ^ . 

Therefore, by the definition of the constant we obtain 

d{Q,T, Lp) ^ diQoF,F-\T),Lp) ■ (detF)^ 
|T|'+p \F^^{T)\^+p ■ {detF)'+p 



\T\ 



Hence, 



dQ,T,p>C;\T\'+-^^ 

miliar 



which completes the proof. □ 

Note that in Section 3 it will be shown that the infimum in ([2]) is achieved on equilateral 
triangles. Thus, to obtain the triangle on which the inequality ^ becomes equality, we should take 
an arbitrary equilateral triangle T, and then the triangle F{T) will be optimal. 

In addition, we shall need the following two lemmas. 

Lemma 5 Let us consider the collection of quadratic forms of type which satisfy the following 
conditions: 

< A < A+, < < 5+, and H = AB-C'^> C+, 
where ^4+, 5+, C+ are some positive numbers. Then for any such form 



Xr^^>\{A+ + B+)-^(^-{A+ + B+)^ -C+>0. 

To prove this lemma observe that the function g{u,v) = u — — v (u > 0, < f < 1) is 
decreasing in u and increasing in v. 

Lemma 6 For the collection of quadratic forms satisfying the assumptions of Lemma \^ the ratio 
of the diameter of the optimal triangle to the square root of the area of this triangle is bounded by 
the constant non-depending on A'^ , and . 

This lemma follows from Lemma 5. 

The next two results will be used in Section 5. 

Lemma 7 Let f G C^{D); H{f] x,y) > > for all {x,y) E D. If n is an arbitrary unit vector 
in the plane, then 



(9) 



d'f 



C 

> := ■ min 



2 lll/xxlloo' 11/^ 



yyiloo 
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Proof. Let n = {u, v) be an arbitrary unit vector in the plane. Then for an arbitrary point 
{x,y) e D 

g^i^' y) = /^^(^' ^)"^ + 2/xy(a;, y)uv + fyy{x, y)v'^. 

Note that functions fxx and fyy have the same sign on D. Without loss of generality we may assume 
that fxx{x, y) > for all (x, y) G D. Since + = 1 then either > ^, or > ^. If > ^ then 

^ ^ 'x,y) = ( ^^^\ '^\ u^ + 2fxy{x,y)uv + fyy{x,y)vA + ( fxx{x,y) - {''^^ '^^ ] u'^ 



9n'^ ' \fyy{x,y) ' ' y V ' fyyi^^y) 

fyyi^^y) J 2||/j^j^||oo 



Similarly, if f ^ > ^ 



^2 i- (J+ 



Thus, we have obtained the desired inequality. □ 

Lemma 8 Let f G C^{D); H{f; x, y) > > for all {x, y) E D. Then for any triangle T 

d{f,T,Lp) >^.d(Q,T,Lp), 

where constant was defined in 

Proof. Let A, B, and C be the vertices of the triangle T. Set 

n{x, y) := /(/, T)(x, y) - f{x, y), (x, y) G T, 

and 

D+ — — 
^{^^ y) '■= ^ ■ {KQ^ T){x, y) - Q{x, y)) , (x, y) G T. 

Obviously, n{x,y) > and m{x,y) > for all {x,y) G T. Moreover, 

n{A) = n{B) = n{C) = m{A) = m{B) = m{C) = 0. 

Let us consider the function 

g{x,y) ■.= n{x,y) -m{x,y), {x,y)eT. 

For the Hessian of g we have 

H{g; X, y) = {-fxx{x, y) + D^) {-fyy{x, y) + D+) - {fly{x, y)) 

= fxx{x, y)fyy{x, y) - fly{x, y) + {D^f - ■ {fxx{x, y) + fyy{x, y)) 

> C+ + {D+f - D+ ■ ifxxix,y) + fyyix,y)) > {D+f. 

Therefore, g is concave on the triangle T, since gxx{x,y) < for all {x,y) G T. It follows that 
g{x,y) > for all {x,y) G T. Thus, n{x,y) > m{x,y) for all {x,y) G T, from which the desired 
inequality easily follows. □ 
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3 The proof of Theorem 2 for 1 < p < oo. 



In this section we shall provide the proof of Theorem 2 in the case when 1 < p < oo and show the 
dependence of the error of interpolation of the quadratic function on a triangle on the geometry of 
the triangle. 

Lemma 9 Let 1 < p < oo. Then for every non- equilateral triangle T there exists a triangle T such 

that _ 

diQ,T,L,) ^ d{Q,T,Lp) 



\l P 



\l p 



Proof. Assume that T is not equilateral, i.e. \AB\ ^ \BC\. Let S be an arbitrary rotation 
of the plane. Then, obviously, Q o S = Q. We may assume that A = (—1, 0), B = (a, b), and 
C = (1, 0), where a e R and 6 > 0. 

Let M = {xm, Um) be the center of the circle circumscribing triangle T, and let R be its radius. 
Clearly, the point M has the following coordinates: 



M= 0, 



+ - 1 
26 



and for the radius R we have: 



Ri' 



2b 



+ 1. 



Then for the error we obtain 

dP{Q,T,Lp) = dP{Qi;--yM),T,Lp) 



(10) 



b^ + a^-l\ ^ 2 



+ - 1 
2b 



dx dy. 



Let T' be the triangle with the vertices A, B' = {—a, b), and C, and let T be the triangle with 
the vertices A, B = (0, b), and C. Obviously, 



(11) 



j^^l i^t'i \t\ 



Let us show that T is the desired triangle. To this end, we consider the linear transformation 
F, which is determined by the matrix 

■l f 
1 

Note that F transforms the triangle T into the triangle T. Hence, we obtain 

dP(Q,T,Lp) = dP(QoF,T,Lp) 



(12) 



1 — u —2--UV 



62 



V — V -\- 



62 + a2 - 1 



dudv. 
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Note that 



d{Q,T,Lp)=d{Q,T',L, 



Therefore, in view of (11) and the triangle inequahty, we obtain 

d(Q,T,Lp) _ d{Q,T,Lp) + d{Q,T',Lp) 



2|T 



> 



1 



1-u' 



62 



V — V + 



du dv 



Since {u,v) E T, we have v < b. Using ( |10[ ) with a = and equahty ( |12[ ), it follows that 

62-1 



d{Q,T, Lp 



1 



IT 



1+ 



> 



IT 



1+ 



1 - - f 2 + 



(iw (if 



rf(Q,T,L,) 



IT 



1+ 



□ 



Lemma 10 Let T be an arbitrary triangle. Then for any < p < oo 



d{Q,T,Lp) > 



diamT 

2^+p/i(T) 



\I p 



where h{T) denotes the minimal height of the triangle T. 

Proof. Let A, B, C he the vertices of the triangle T, and let \AB\ = diamT. We may assume 
that A = (— a,0), B = (a,0) and C = {b,h{T)), where 2a = \AB\ and b G {—a, a). Consequently, 
a ■ h{T) = 1. Without loss of generality we may assume that 6 > 0. Let us consider the trapezium 
PQST with vertices P = (^, ^), Q = ^), S = (|,0), and T = (-|,0). Note that the 

area of the trapezium PQRS is equal to ||T|. 

Let M = {xm^Um) be the center of the circumscribed circle of the triangle T, and let R be its 
radius. Obviously, 



M 



and i?^ = + y^j. 



Then 



Recall that 



L := min{i?2 - \MS\^; W - |MTn W - |MP|^; W - \MQ\^} > 



16' 



dP{Q,T,Lp)= I - (x - xmY - (y - VAiff dxdy. 



Moreover, 

dP{Q,T,Lp) > 



\Fi^ — {x — xmY — {y — yM)^'f dx dy > U' I dxdy> 



fl^p . |r| 

2 ■ 16P 



PQST 



PQST 
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Hence, 



dP{Q,T,L,)> 



\AB\^P\T\ \AB\P\T\P+^ 



26p^ 



□ 



25p+ihP{T) 

In view of Lemma 10 we can derive that the infimum in the right hand side of ^ exists. 
Moreover, by Lemma [Opthis infimum is attained on the equilateral triangles only. 
The proof of Theorem 2 in the case p G (0, 1) can be found in the Appendix. 



4 Error of interpolation of functions by linear splines: 
upper estimate. 

In this section we shall show that 



limsup N ■ Rn if, Lp^n) < ^ ( I H{f;x,y)^Q{x,y)^^dxdy 



p+i 

V 



In order to do so we are going to construct a sequence of triangulations {A^}^^^ such that 



hmsupiV- ||/-s(/,A^)|U^^ [ H{f-x,y)^Q{x,y)^^dxdy) 

Af-5-oo \J D J 

For a fixed e G (0,1) and for every G N we define 



£+1 



(13) niN = TaM{e) := min <^ m > : — ^ w ( — ) < 



m? \m I ~ N 



where w(5) is the function defined in (|4]). 

Observe that m^v— j-csoasA^— j-oo. In addition, note that 

A^ 

(14) —2 )• oo, A^ -> cx), 



i.e. niN = o (^V^) as A^ — cxd (see, for instance. Section 4 from [1]). 

Let us divide the square D into squares with the side length equal to ^ that have sides parallel 
to the sides of D, and denote the resulting squares by := Df{e), i = 1, . . . , mif, enumerated in 
an arbitrary order. 

Let {xf, yf), i = 1, . . . , m^, be the center of the square Df. In addition, for every i = 1, . . . , 

set 
and 

:= sup il{x,y). 

Note that 

if(xf , y^) := H{f; xf , ) = 4(Af i?f - (Cf )^), ^ = l,...,m%. 
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Set 

Qf (x, y) := Afx' + 2C^ xy + 5,V, (x, v) e M'- 

Now for the fixed e and for all N large enough we will construct an appropriate triangulation 
/\n{^) of D consisting of triangles. To this end, we shall construct the triangulations : = 
A^(e) of squares Df depending on the eigenvalues of the quadratic forms Qf , i = 1, . . . ,m%. 
After this we shall "glue" these triangulations to obtain the triangulation A^is) of D. 

Everywhere below fci, ^2, . . . , stand for constants independent of and e. 

Note that for an arbitrary triangle T E D and for an arbitrary continuous function g G C{D) 

(15) d{g,T,Lp^n) < \\n\\^\\g\\L^^T)\T\p . 

For every i = 1, . . . , let 

fN,i{x^ y) = P^if; X, y; xf , ), 

where P2{f', x, y; a, b) is the Taylor polynomial of degree 2 of / constructed at the point (a, b). 
For every i = 1, . . . , mj^ set 



= (.) 



p j\r p 



+1. 



Observe that all nf^ — )■ oo when N ^ oo. This follows from the obvious estimate 



(16) 



min y)}p^ 

(x,y)eD 



oo 



1, 



combined with (14), and min VL[x,y) > 0. 

I (a;,y)GD ^ 

For an arbitrary triangle T C let us consider the triangle T which is symmetric to T with 
respect to the midpoint of one of the sides of T. Therefore, 11 := T U T is a parallelogram. By 
tilling the plane with the help of 11 we shall obtain a triangulation of the plane G{T). 



Given nf for each square Dj 
as follows: 



N 



,mj^, we construct the triangulations Ajy of the square 



1. We consider transformations Fi and F2 of the form ([T]) and ^ respectively, corresponding to 
the quadratic function Qf{x,y). Set F := Fio F2. 

2. We take an arbitrary triangle T which solves problem ([2]), and consider the triangle F{T). 

3. Next we define Tl^ to be a re-scaling of F{T) so that 

1 



N\ 



2 N ' 
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4. With the help of the triangle Tj^ we generate the triangulation GiTj^) of the plane as described 
above. 

5. Triangles from GiTj^) which lie completely in the square Df we include into triangulation 
A^. Let us consider the triangles from G{T^) having common points with the boundary 
of Df . The intersection of an arbitrary such triangle with is a polygon with at most 7 
vertices. We triangulate this polygon into at most 5 triangles without adding new vertices. 
All triangles which we have obtained are included into triangulation A^. 

Now let us "glue" triangulations A^y. We shall do it according to the following algorithm. 

1. For every i = 1, . . . , m% let Wf^ be the set of the vertices of triangulation AJ^ which lie on 
the boundary of Df^ . 

2. For arbitrary i, j = 1, . . . , m%, i ^ j, set Sij = Df fl Dj^. 

3. Let Sij 7^ 0. For every triangle T G A'^ if T fl Sij is a non-empty segment then we subdivide 
this triangle by joining the vertices of T and Wj^ (1 T. Then we include all these triangles in 
the triangulation A]\[{e). 

Therefore, for every e G (0, 1) there exists N{e) G N such that for all N > N{e) we present a 
construction of the triangulation A^ie)- 

Obviously, for every i = 1, . . . ,m% the conditions of Lemma 6 are satisfied for the quadratic 
form Qf. Thus there exists a constant ci, independent of N, such that for every i = 1, . . . , 171% and 
for any triangle T from the triangulation A^y we have 



N 



(17) diamT < cn 

for all large enough. 

For all i = l,...,m^ denote by the number of triangles from An{b) which lie in Df. 
In addition, we denote by A^i = Ni{e) the number of triangles in the triangulation A]\f{e), and 
let A^i = Ni{e) be the number of triangles from the triangulation Aj\f{e) which have nonempty 



m 



intersection with |J dDf . 

i=l 



From (17) it follows that for all large enough (without loss of generality we may assume that 



this is true for all A^ > N{e)) 

(18) nf<(l + £)nf for all i = l,...,m%, 

A^i < (1 - £)N and Ni < eN. 

In view of construction of the triangulation Ai^{e), for every triangle T G Ai^j{e) which lies 
completely in int Df (the interior of the set Df) we have 



(19) d^{fN,i, T, L,) = <FUN,i, T^, L,) = (^^^ 



13 



Then for the error of interpolating fj^^i on the square Df at the vertices of A^y we obtain that 

TCint 



Hence, in view of (18) and (19) 



'"AT /(7+ \ ^ ™^ 



rN 



n. 



N 



Using the definition of nf we obtain 



Vi < 



2 I -^2p+2 



N i=l 



1 + e 



i=i 

P AT p 

N{1 - e)H{xf ,yf)^(^^ )^ 



P+i 



1 + e 



P+i 



(1 - e)PiVP I m' 



1 

% i=i 



Hence, by Riemann integrabihty of functions H{f] x,y) and fl{x,y) 



1 i\ 

lim y/7(xf,yf)^(nf)^*T = limy \Df\H{x' 



I 2(p+l) 



|p+i 



D 



Hence, for all large enough we obtain 



(1 + he) (^j^ H{f- X, y) y) dxdy^ 



p+i 



Note that for all e e (0, 1) and for all > N{e) we have 

(20) RN{f,Lp^n) < ||/-s(/,A^(£))||p,o. 



Let us consider the case 1 < p < oo. Observe that 



(21) 



ii/-.(/,A^(e))ii^,^ = J2\\f - s{f,A^em;,n<J2\\\f - 



i=l 



i=l 



+ \\fN,i - s{fN,h AN{e))\\p,n + \\s{fN,h A7v(£:)) - AN{e))\\p^n 
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Obviously, for every i = 1, . . . ,m 



N 



|s(/Ar,i, Aiv(£)) - s{f,ANi£))\\p,n < -^\\f - fN,i\\oo,n- 



m 



N 



and 



11/ — fN,i\\p,n < — n-\\f — /Af,i||oo,C- 



m 



N 



In addition, by Lemma [T] and the definition of m^v we liave 



2||^^||oo 

11/ — /iV,j||oo,n < 



m 



N 



-) < >\ 



Thus, 



Hence, 



'"JV 

ll/-^(/,A^(5))||^,^<5^[||/;v.-3(/;v,,,,A;v(e))|U: 



i=l 



2e 



||/-.(/,A^(5) 



\p,n 



< 



i=l 



< (1 + he) J2 Un, - s{fN,, A^(£))||^^^ =: (1 + he)V,. 



i=l 



p+1 



Therefore, there exists N{e) > N{e) such that for all > N{e) 

l + kis) i^j H{f;x,y)^(^il{x,y)^dxdy^ 
Note that in the case p G (0, 1) instead of (21) we should use the following implication of ([T]): 



(22) ||/-.(/,A^(e))||^,^< (^'' 



||/-^(/,A^(e))||^,^ 



p 



i=l 



i=l 



+ \\fN, - s{fN„ A^(£))||^_^ + \\sifN„ A^(£)) - S{f, A^(£ 



\p,n- 



The rest of the proof is analogous to the case 1 < p < oo. 

Let 1 > £1 > £2 > • • • be a decreasing sequence of positive numbers which tends to zero as 
/c — )■ oo. Without loss of generality we may assume that {N{ek)}'^=i increasing sequence. 

Then set 

A^:=A^(£,) if Ni6k)<N<N{ek+i), keN, 
and let A^ be an arbitrary triangulation of D if 1 < < N[ei). 



Therefore, in view of the definition of -RAr(/, Lp Q) and inequalities (20) and (22), for all < p < 
oo and for every N{ek) < N < N{ek+i) we have 



RN{f,Lp,n) < \\f- s{f, A%)\\,^n < ^(1 + to) (^j^H{f-x,y)^n{x,y)^^dxdy 
Since — )■ 0, as /c — )■ oo, we obtain the desired estimate. □ 



£+1 
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5 Error of interpolation of functions by linear splines: 
lower estimate. 

For an arbitrary triangle T in the plane denote by h{T), diamT, U{T) and |T| the minimal height 
of T, the length of the longest side of T, an arbitrary vertex of the longest side of T, and the area 
of T, respectively. 

To prove the lower estimate we need to show that 



(23) \imMNRN{f,Lp,n)>^{ I H{f;x,y)^n{x,y)^^dxdy 



D 



E+l 

V 



Let VL = min Q{x,y), and let {Ajvj^^j^ be an arbitrary sequence of triangulations of D. 

ix,y)eD 

For any e > let us consider the following sets: 

Me) :={.€{!,..., iV} : < -^^(diamTn), 

where 

8||f^||oo 

c : = 



^Ov/WTfW)) 

and 



Jiv(e) := 



{^e{l,...,iV} : diam7;^>^}. 

Note that the sets In{£) and Jn^^) niight have a nonempty intersection. 

First, assume that there exists > and a subsequence {Nk}'^i of positive integers such that 

(24) \Tl'''\>^o- 
Let us show that in this case liminf A^fc||/ — s{f, Ajv^)||p^Q = oo. 

Indeed, in view of Lemma [t] there exists a constant such that the derivative of the function 
/ in an arbitrary unit direction is bounded away from zero by D^. Thus, by Lemma [s] we obtain 
that 

\\f-s{f,/^NM;,n > E d%f,T!''',L,,n) 

iei'jv^(eo)\^]Vt.(eo) 

ie/iVj. (eo)\-/jVj.(eo) 

Hence, by Lemma 10 and the definition of iNki^o) we have 

'QD+V 1 ^ /diamr^'= 



> 



\ 2c J 25p+i ^ wPfdiamT.^M' 
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Let card (F) denote the number of elements in the finite set F. Since diam (T^ ) < applying 
the Jensen inequality and (24), we obtain that 

Eonn 



\\f-s{f,Ar,,)\\l^ 



> 



2c ■ LO 




25p+l 



E 



ie-r]Vj.(eo)\^JVj.(£o) 



\ i I 

p+1 



> 



> 




E 

^&lN^{eo)\JNf, (eo) 



card(/w,(5o) \ JN,{eo)) 



The last inequality implies that liminf A^a;||/ — A7Vfe)||pQ = oo. 

Assume now that there exists a number Eq > and a subsequence {iV^}^^ of positive integers 
such that 

«eJjVj.(£o) 

As in the previous case we are going to show that liminf Nk\\f — s{f, AN^)\\p^n = oo. 

fc— >oo ' 



Similarly to above we can derive that 

/'nD+\^ 1 



ieJjv,, (so) 



diam T- 



Therefore, in view of |T| = |diamT/i(T) for an arbitrary triangle T and the definition of Jni,{£o), 
we have 

/ iGJivJeo) 



> 



> 



This implies that liminf A^^ ■ ||/ — A 



fc— >oo 



Nk)\\pn 



oo. 



Hence, in order to find a sequence of triangulations {AatI^^^i which provides the lowest value 
of liminf A^ ■ ||/ — Ajv)||pf7, we should investigate only the sequences of triangulations with the 

TV— >-oo ' 

following property: for an arbitrary £ > there exists N{e) such that for all A^ > N{e) 

ie/jv(e)UJiv(£) 

Let {Aat}^^]^ be a sequence of triangulations satisfying this condition. For every z = 1, . . . , A^ let 
/jv,i be the second degree Taylor polynomial of / at the point U{Tl^). Then in view of Lemma 1 
we obtain that 

11/ - fN,\\L^(Tn < 2(diam7;^)2a;(diam7;^) 
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for alH = 1, . . . , A^. Therefore, 

ll/-/iv.llL„„(i;-) < ||f^||oo||/-/Ar,|k^(T^)|7;^'|^ <2||0||oo(diain7;^^)^a;(diam7;^^)|^ 

and 

ll^a,?;"^) < 2||n|U(diam7;^)2a;(diam7;^)|7;^|^. 

Consequently, for all i ^ In{^) 

11/ - /iV,i||Lj,,n('Z;^) + - ^(/iV,i>^^)llLj,,o(Tf) 

< 4||Q|U(diam7;^)2a;(diam7;^)|7;^|^ 

< 4||0|U(diam7;")^|7;"|^ ^^;^|^^^^^, ^ ^^JPU^^n^..^ 



< e-^-^l^H{f-U{Tn).\Tl^rl. 
Set Of := inf Vtix, y), i = 1, . . . , N . Obviously, 

N 

(25) \\f - s{f,Ar,)\\l^^J2d^if,T!',L,,n) > ^ <F{f,T^,L,,a). 

i=l i07jv(£)UJjv(£) 

Applying the triangle inequality and the defintion of sets In{s:) and Jn{£^), for all 1 < p < oo 
and i — 1, . . . , N we have 

dP{f,Tf,L,,a) > {d{fM,i,T^,L,,a) 

- 11/ - /iV,i||Lp,n(Tf ) - IIH/>^^) - H/iV,j,7^^)||Lp,n(Tf 

- 4||l]|U(diam7;^)2a;(diam7;^)|7;^|^^)'' 

- ((1 - ^) • ^(^'')) • Of ■ 

Similarly, for p e (0, 1) and i = 1, . . . , A^" we have 

i;^, L,,^) > (^(1 - he) ■ {^jHUmTf)) ■ 9.^y\Tn^^\ 



p 
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Therefore, in view of (25), for all p G (0, oo) 



> 



i0N{e)UJN{s) 

(1 - k,e)C^ 



2{N-\lN{e)uMe)\) 



p+i 



X 



> 



j^/jv(e)UJiv(e) 

{l-kre)C;Y 



Ni 



p+1 



2N 



Let us divide each triangle , i G /tvI^) U JNi^), into nf^ triangles T^^, j = 1, . . . , raf^, enumerated 
in an arbitrary order, such that diam (T.^) — )■ as — )■ oo for all j = 1, . . . , nf^ . For every 
i ^ /Ar(£) U JAr(e) Set nf' = 1 and T^^ = Tf. In addition, set Oi^- := inf il{x,y) for all 



N nf 



j = 1, . . . ,nf^ and z = 1, . . . , A^. Note, that [J IJ Tfj = and for every i and j we have that 



i=i j=i 



diam (Tfj) — )■ as — )■ oo. Then 



ll/-^(/,A^)||^,^ 



> 



{l-kje)C+\ 



2N 



P N K 

^ * i=i i=i 



2(P+1) if^M ~^-JL^\rpN \ 



P+1 



(ll^lloovW 



I p+1 



ie/iv(e)UJjv(e) 



> 



p+1 



2(1 + k^e)N 



Therefore, 



li,„mfiV||/-.(/.A.)|U„>^l^ 



and since e is arbitrary we obtain the desired inequality. □ 



6 Appendix 



In this section we shall show that for < p < 1 



inf 

T 



d{Q,T,L,) 



4 



3^3/ 



TT 



P+1 



1 

6 J x(l 

1/2 



x^)^ arccos — dx 
' 2x 
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Note that 

(26) C; = ^iM^t/(Q,T,L,). 

By M and R we shall denote the center and radius of the circle circumscribing triangle T, respec- 
tively. We may assume that the point M is at the origin. Then 




Figure 1: The triangles Ta = ABMC, Tb = ACM A, and Tc = AAMB, in the cases of acute 
(left) and obtuse (right) triangle 

Let A, B, C be the vertices of the triangle T. By the same letters we shall denote the angles 
corresponding to the vertices of the triangle T. Without loss of generality, we may assume that 
\AB\ = diamT, and, consequently, C > A and C > B. Let T^, Tb, and Tc be triangles obtained 
by joining the point M with the vertices A, B and C (see Figure [T]). 

Set 

m{A) := [ [{R^ -x^ -y^ydxdy, 



Ta 



m{B) := jj{R^-x^-y^fdxdy, 

Tb 

m{C) := [[ {R^ - - y^ dx dy. 



Tc 

Note that in the case of non-obtuse triangles 

dP{Q, T, Lp) = m{A) + m{B) + m(C), 
and for obtuse triangles (C > |) we have 

dP(Q, T, Lp) = m{A) + m{B) - m{iT - C). 
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Let the triangle T\ be homothetic to the triangle T4 with the side, corresponding to the side OB, 
equal to 1. From the definition of m{A) it follows that 

(27) m{A) = i?2P+2 y2y ^ 2i?2p+2 _ i^A)] , 



2p + 2 



where 

1 



1{A) = j x(l — x^)^ arccos ^^^^^^^ — ^dx. 

cos A 



Note that /(O) = 0, /(7r/2) = and 

=7(p)(sinA)2f+2, 

where 7(p) = [p+ 1, 




Figure 2: The domains Fi (left) and F2 (right) 

Every triangle with unit area is uniquely determined by the pair of its angles {A,B), which 
can be considered as the point in the plane with coordinates {A,B). Let Fi and F2 be the sets of 
points, which corresponds to the sets of all acute and all obtuse triangles respectively (see Figure [2]). 
Obviously, Fi and F2 are open sets in the usual topology of M^. 

In view of Lemma 10, the infimum in (26) is achieved at some triangle Ti. Let Ai, Bi and Cibe 
the angles of the triangle Ti. 

First, assume that the triangle Ti is acute. Taking into account the expression (27), in this case 
we have 



dP{Q, T, Lp) = m{A) + m{B) + m{C) = 2R^p 



+2 



TC 



2p + 2 



1{A) - 1{B) - 1{C) 



Since ITI = 1, we obtain 



sin 2A + sin 2B + sin 2C 
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Hence, 



(28) 



Z{A,B) 



2P+2 



2p+2 



1{A) - 1{B) - 1{C) 



(sin 2A + sin 2B + sin 2C)p+i ' 



Therefore, by the choice of the triangle Ti, 

^(Ai,5i; 



inf ZiA.B). 

{A,B)eFi 



Since Ti is an acute triangle, the necessary conditions of extremum have to be satisfied: 



dZ_ 
dA 



V{C) - V{A) 



(29) 



and 



(30) 



(sin 2A + sin 2B + sin 2C)p+i 
(2p + 2)(cos 2A-COS 2C) 



2p+2 



1{A) - 1{B) - 1{C) 



(sin 2A + sin 2B + sin 2C)p+2 



9Z 



/'(C) - l'{B) 



(sin 2A + sin 2B + sin 2C)p+i 
(2p + 2) (cos 25 -cos 2C) 



2p+2 



1{A) - 1{B) - 1{C) 



(sin 2A + sin 25 + sin 2C)p+'^ 







From equations (|29|) - (|30|) we derive that either the angle Ci is equal to one of the other angles of 

/'(Ci) - l\B,) 



the triangle Ti (Ci = or Ci = -Bi) or 

l\C,) - l\A,) 



(sinCi)2 - (sin/li)2 (sinCi)2 - (sin5i)2" 
From the last equation with the help of the expression for l'{A) we have: 

(sin Ci)2p+2 - (sin Ai)2p+2 (sin CifP+^ - (sin BifP+^ 



(sinC/ 



(sin Ai) 



(sinCi)2- (sin5i)2 



which can be true only if Ai = Bi. Therefore, if Ti is acute then it is isosceles. 

Now, if we assume that Ti is obtuse then we similarly obtain that Ti should be isosceles. 
Therefore, the triangle Ti solving (26) must be either right or isosceles. 

Now let us consider the case of right triangles. Then C = |. We shall show that the triangle 
solving the problem 



c/(Q, T, Lp) ^ inf , |r| 



1, C 



2' 



is isosceles. Indeed, in this case B 



dP{Q,T,Lp) = m{A)+m 



A and 



TT 



TT 



A 



m 



TT 



2 ■ 



4p + 4 



(sin 2A) 



P+i 
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Lemma 11 The function 



(31) 



4p + 



(sm2A)p+i 
Z5 non-increasing on the interval (0, 7r/4]. 

Proof. Let us consider the derivative of the function L: 



L'(A) = 



[/'(7r/2 - A) - /'(A)] sin 2 A - {2p + 2) cos 2 A 



Ap + A 



1{A) - l{7i/2 - A) 



(sin 2^) 



p+2 



Taking into account the expression for I' (A), in order to show that L'{A) < on the interval (0, n/A) 
it suffices to prove that 



r{A) := 7(p)[(cos/l)2f+2 _ ^smAfP+^] tan2yl - (2p + 2) 

on the interval A e (0, 7r/4). 

Let us consider the derivative r'{A): 



71 



Ap + A 



- 1{A) - /(7r/2 - A) 



< 



"'^"^^ ^ 7(p){(^^^[(cosA)2^+2-(sinA)2^+2] 

- {2p + 2)tan 2A sin A cos A[(sin A)^^ + (cos A)^p] 

- {2p + 2)[(cos^)2f+2 - (sin^)2f+2]|, 



where 'y{p) = {p+ 1, |). 

Set t — tan A, and then t e (0, 1]. Obviously, 



r'{A) = 



l{p)z {t) 
(cos A) 



2p+2' 



where 



z{t) = peP+^ - 2t2p+4 _ (p + 2)t2p+2 + (p + 2)t^ + 2t^ - p. 

To prove that r'{A) < for all A e (O, |] we shall need the following proposition. 
Proposition 1. Let p E (0, 1). Then the function 

z{t) = pt^P^"" - 2t^P+^ -{p + 2)t^P+^ + {p + 2)i^ + 2t^-p 



TT 



is non-positive for all t e [0, 1]. 

Thus, we obtain that r'{A) < for all A e l^*^' ^ • Consequently, the function r{A) is non- 
increasing, and hence r{A) < r(+0) = 0. Hence, we have proved that the function L(A) is non- 
increasing on the interval (0,7r/4], i.e. 



L{A) > L{n/A) 
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for all A G (0,7r/4). □ 

Note that it can be shown similarly that the function L{A) (see (31)) is non-decreasing on the 

interval — , — V Therefore, 
4' 2/ 



inf d{Q, T, Lp) = d{Q, T2, Lp), 
|r|=i,c=f 

where T2 is the isosceles right triangle. 

Hence, we have proved that the triangle Ti must be isosceles. Without loss of generality we may 
assume that A = B. Let us show that Ti is equilateral. To this end let us consider the following 
problem 



(32) 



d{Q, T, Lp) inf, |T| = 1, A = B. 



Let us consider now the case when the triangle T is isosceles non-obtuse triangle. Then A G 
|,0,and 

21{A) - 1(71 - 2A) 



dP(Q, T, Lp) = 2m(A) + mi-n - 2A) = 2^+^ ■ ^ — — 

' V ;^ V ; (2sin2A-sin4A)p+i 



Lemma 12 The function 



SiA) 



2p+2 



- 21{A) - l{u - 2A) 



(2sin2A-sin4A)p+i 
is non-increasing for A G [7r/4, 7r/3] and is non- decreasing for A G [7r/3, 7r/2). 



Proof. Let us consider the derivative of the function 5* 

8cosA(sinA)3[2/'(7r - 2A) - 21' {A)] 



S'{A) 



(33) 



(8cosA(sinA)3)p+2 
8(p+ l)sinAsin3A 



2p+2 



21{A) - /(tt - 2A) 



(8cosA(sin A)3)p+2 



Obviously, the denominator of the right hand side of (33) is positive in the considered region. In 
what follows we shall consider only the numerator of (33). Then showing that S'{A) < 0, when 
A G [7r/4,7r/3] and S'{A) > when A G [7r/3,7r/2), is equivalent to proving the inequality: 



q{A) := 27(p) 



cos A sin A 
3-4(sinA)2 



;sin2A) 



2p+2 



(sin A) 



(p+1) 



TT 



2p + 2 



21{A) - l{n - 2A) 



2p+2 



< 



for all A G [77/4,77/2). 
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It is easy to see that g'(7r/2 — 0) = and 

27(p) 



q\A) = 



[3-4(sinA)2]2 



{ [(sin2A)2f+2-(sinA)2f+2] 



X 



cos2A[?, - 4(sinA)2j + S{cosAf{smAf 
- (p + l)(^3-4(sinA)2)^ 
+ 2(p + 1) cos A sin A 3 - 4(sin Af 

X [2 cos 2^(sin 2yl)2p+i - cos ^(sin AfP+^] | . 



Obviously, 



where 



l{p) 



[3-4(sin>l)^ 



i = 2cosA e (0,^2], 



Z{t) := (i2p+2_l). (^2_2)(^2_;L)+i2(4-i2)_2(p+l)(t2_l)2 

+ (p + l)i(i2 - 1) [2(^2 _ 2)i2P+l _ i 

= -{2p + l)t2p+4 + 2(p + 2)i2p+2 + (p + i)t4 _ (3p + 4)t2 ^ 2p. 

We shall need the following proposition to show that z{t) changes its sign exactly once. 
Proposition 2. Let p e (0, 1). Then the function 

z{t) = -(2p + + 2(p + 2)t2f+2 + (p + 1)^4 _ + 4)^2 ^ 2p 

has exactly one point of sign change (from positive to negative) on the segment [0, 2], and this point 
is located inside the interval (0, 1) . 

Therefore, z{t) changes its sign on the segment [0, \/2] exactly once from positive to negative. 
Since variable A increases when variable t decreases, the function q'{A) changes sign exactly once 
as well, from negative to positive. Hence, 

q{A) < min{g(7r/4); q{n/2 - 0)} = 0, 

whenever g(7r/4) = 7(p)(l — 2~'p~^) — (p + 1) — 2/(7r/4) < (this will be shown in the proof 

of Lemma 13). Hence, we have proved that S'{A) < for ^ e [7r/4,7r/3], and S'{A) > for 
A e [7r/3,7r/2). □ 

Let us turn now to the case of isosceles obtuse triangles T. In this case A e (O, |), and 



(F(g, T, Lp) = 2m{A) - m{2A) = 2^+^ 
Lemma 13 The function 



-21{A) + 1{2A) 



S{A) 

is non-increasing on the interval (O, |] . 



(2 sin 2A - sin AA)p+^' 
2l{A) + l{2A) 



(2 sin 2 A - sin 4A)p+i 
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Proof. Let us consider the derivative of the function 5* 

8cosA(sinA)3[2/'(2A) - 21' (A)] 



S'iA) 



(8cosA(sinA)3)P+2 
8(p + 1) sin A sin 3A[-2l{A) + 1{2A)] 
(8cosA(sin Ay)P+^ ' 



We shall show that S'{A) < when A G (0, tt/A] or, equivalently, that q{A) < for all A G (0, 7r/4] 
where 

... , , cos A sin A 

q{A) := 27(p), 



(sin2A) 
3-4(sinA)2L' ^ 

- {p+l)[-2l{A) + l{2A)]. 
The derivative of q{A) can be written as: 

27(p) 



2p+2 



(sin A) 



2p+2 



[3-4(sinA)2]' 



I [(sin2A)2p+2_(sinA)2p+2] 



cos2A(^3 - 4(sinA)2j + 8(cos A)2(sin A)^ 



7(p) 



- (p+ 1)(^3 -4(sinA)2^ 
+ 2(p + 1) cos A sin A 3 - 4(sin A)^ 

X [2 cos 2A(sin 2A)2p+i - cos A(sin A)^p+'^] | 



[3-4(sinA)2] 
where t = 2 cos A G [72, 2) and 

z{t) = -{2p + + 2{p + 2)t2p+2 + + i)t4 _ + 4)^2 ^ 

In view of Proposition 2, it follows that z{t) < for all t G [-\/2, 2). Hence, q'{A) < on the 
segment (0, vr/4]. Therefore, q{A) < g(+0) = (in particular, g(vr/4) = g(vr/4) < 0) and S{A) is 
non-increasing on the segment (0, 7r/4]. □ 

Combining Lemmas 12 and 13 we conclude that 

inf d{Q, T, Lp) = d{Q, T3, Lp), 



where T3 is equilateral triangle. Moreover, T3 is the only triangle which solves (32). Since the 



triangle Ti solves (26) and is isosceles, we obtain that Ti = T3, which finishes the proof. 
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